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s ! ■ We show that the combined action of diffraction and convection (walk-off) in wave mixing pro- 

j_( ■ cesses leads to a nonlinear-symmetry-breaking in the generated traveling waves. The dynamics 

near to threshold is reduced to a Ginzburg-Landau model, showing an original dependence of the 
, nonlinear self-coupling term on the convection. Analytical expressions of the intensity and the ve- 

locity of traveling waves emphasize the utmost importance of convection in this phenomenon. These 
predictions are in excellent agreement with the numerical solutions of the full dynamical model. 

m : 

PACS numbers: 89.75.Kd,47.54.+r,42.65 Sf 

Qh \ There is currently a considerable interest in understandingthe role of convection in pattern forming systems 

in such diverse fields as hydrodynamics Q , plasma physics p , traffic flow Q , crystal growth Q , and nonlinear optics 
Q . The most important and common result, in these studies, is that the pattern selection in spatially extended systems 
is dramatically affected by the breaking of the reflection symmetry (r — ► —r) due to the presence of convection. This 
is a linear symmetry breaking where convection terms are generally considered to have the only effect to induce 
a traveling character to selected patterns, and to induce a peculiar regime of convective instability. It has been 
largely studied how the existence, the type and the dynamics of the pattern selected are closely related to the linear 
\ transition from convective instability (propagation overcomes amplification of perturbations) to absolute instability 
\Q i (amplification dominates). 

In contrast, in this paper, we discuss an unforeseen effect of convection in the dynamics of spatially extended 
systems that does not rely on such a transition. Here, we show how convection, that is a linear phenomenon, actually 
modifies the intrinsic nonlinearity of the system. More precisely, we show that convection affects the nonlinear 
modes interaction at onset of the instability leading to a nonlinear symmetry breaking in the generated, otherwise 
symmetrical, traveling waves. Our analytical description of this mechanism, based on the amplitude equation of 
the degenerate optical parametric oscillator, demonstrates an original dependence of the nonlinear self-coupling term 
upon convection. We consider a degenerate optical parametric oscillator (DOPO) because convection (walk-off) arises 
naturally from the birefringence of the crystal that composes this device. Moreover these devices are at the basis of 
^ . interesting quantum phenomena, stemming from their quadratic non-linearity, as for instance entanglement between 
off- axis modes 0. It has been shown that the walk-off strongly influences such twin beams correlations: in the 
convective regime (induced by walk-off) the entanglement is destroyed by macroscopic amplification of quantum noise 
|lfjj . Increasing the pump intensity, an absolutely stable traveling pattern arises in the signal and the entanglement 
is restored. Still, important walk-off effects are observed, as one of the twin beams is more intense and it fluctuates 
more than the opposite one Although, we present our investigations in the context of optics, we believe that our 
result is generic for spatially extended systems with convection and characterizes the key role of convection in the 
nonlinear dynamics of such systems. 

We start from the description of a type I phase-matched DOPO in the mean-field approximation including diffraction 
and walk-off 01 

d t A p = j p [-(l + iA p )A p + ia p V 2 ± A p -A 2 s + E ] (1) 
d t A s = y s [-(l + iA s )A s +ia s V 2 L A s + A p A* s -a s d x A s \ 

where A p and A s are the normalized slowly varying envelopes for pump and signal fields, respectively. The parameters 
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A p . s , 7p,s, and a PtS are the detunings, the cavity decay rates and the diffraction coefficients, rcspectiveiy. Eq is the 
normalized external pump and a s is the signal walk-off coefficient that characterizes convection in this system. Wc 
stress that the walk-off cannot be eliminated by a change of reference frame, being relative between pump and signal. 

Both convective and absolute instabilities have been reported for the stationary solution A p = Eq/(1 + iA p ), A s = 
[l2^ . We just recall that for the case A s < and A s A p — 1 < 0, which we are interested in, degenerate OPOs exhibit a 

supercritical bifurcation at Eq = + A 2 . This is the linear threshold at which stationary homogeneous solutions 

become unstable to traveling wave perturbations with wave vectors k = (k x ,k y ) with modulus k = k c = y/—A s /a s 
and frequency tu c = —^ s a s k c . Under periodic boundary conditions the convective instability is suppressed and the 
traveling rolls arising at the signal generation threshold are absolutely stable. 

To study the role of convection in the nonlinear symmetry breaking of the generated traveling waves, and to keep 
mathematics as simple as possible, we perform the reduction of the model Q into a single Ginzburg-Landau (GL) 
equation valid close to threshold of the DOPO emission. In the sequel we set fi p = Eq/(1 + iA p ) with /x = \fi p \ = 

Eo/Jl + A|, and B = A p - \x c where \x c = 1. We expand the signal as A s = eA^ + e 2 A^ + e 3 A (3) + ... with a 

similar expansion of the pump. The small parameter e measures the distance from the DOPO emission threshold: 
e 2 = fi — /i c . Setting A^ = A exp i(uj c To + k c Xo) + A* exp — i(u) c Tq + k c Xo) and applying the solvability condition at 
order 0(e 3 ) we get the following amplitude equation, which describes the evolution of the signal written in the scaled 
time r = ^ s t, and in the variable S = eA 

d T S + a s d x S = (/x - 1)5 - 2a s A s C 2 xy S - (3 \S\ 2 S (2) 

with 

2 C 2 + 1 + D 2 + iD(C 2 - 1 - D 2 ) 
/_ 1 + A2 + [1 + (C + D) 2 \[l + (C - D) 2 ] { ) 

We have set C xy — [d x +d 2 /2ik c ) in Eq. J2J and C = A p + Aa p k 2 and D = 2uj c /j p = —2a s (j s /-f p )k c in the expression 
of the nonlinear self-coupling coefficient (3. In absence of convection (a s = 0), the parameter D vanishes and the 
expression of /3 greatly simplifies to j3 = 2/(1 + A^) + 1/(1 + C 2 ) and stationary rolls arise in the signal profile [T^ . 
We show here that the presence of convection drastically affects the pattern formation mechanism with respect to 
both linear and, more importantly, nonlinear dynamics. Indeed, the most important result is that the parameter 
D, that characterizes convection, strongly modifies the nonlinear self-coupling term (3. It affects the saturation term 
[Re(/3)] and induces intrinsic nonlinear phase modulations [Im(/3)]. This result is in contrast with almost all previous 
studies of model equations describing the near-threshold dynamics such as the Ginzburg-Landau or Swift-Hohcnbcrg 
equations where the convection only yields to the propagation term of Eq. J2J. No study has yet been reported, to 
our best knowledge, on any dependence of the nonlinear coefficient j3 upon convection (a s ). 

At this stage one has to notice that the presence of convection via the non zero Im(/3) breaks the well-known 
variational form of the GL Eq. @ in its one-dimensional version (where C xy = d x ), since all the remaining coefficients 
are real 01- A s a consequence Eq. obviously, cannot exhibit stationary homogeneous solutions (or stationary 
rolls). Moreover, the non variational effect leads to (i) a symmetry breaking in the opposite traveling waves and 
(ii) an excess velocity (with respect to convection velocity) in these waves stemming from the nonlinear frequency 
modulation. Both points are analytically characterized below. 

Let us find the solutions of Eq. J2J, corresponding to the nonlinear saturated selected modes, in the form S s t = 
So exp i(flr + kx) . They read |5 | = (m~ 1 + 2a s A s k 2 )/ Re(/3) and Q = — a s k — Im(/3)|5o| 2 , and represent the leading 
contribution to the fundamental modes (±/c c ). This leading contribution is not sufficient since the total intensity of 
each mode is now fixed during their nonlinear interaction induced by the convection. For our purpose to characterize 
the nonlinear symmetry breaking, we need to take into account the contributions up to the third order in e. This can 
be achieved by solving the hierarchy of the inhomogencous linear problems, at each order in e, by means of Frcdholm 
alternative. After lengthy but straightforward calculations, we get the solution 

A s = [1 + (l/2)F 3 \S st \ 2 ]S st expi(uj c t + k c x) 

+ [l-(l/2)F 3 *|S st |XW-*(^ + M (4) 

where F 3 is defined as Re(F 3 ) = 2CD/Den and Im(F 3 ) = 2A p /(l + A 2 ) + C(l + C 2 - D 2 )/Den with Den = 
[1 + (C + D) 2 }[1 + (C — D) 2 ]. Note that the spatial modulations of these traveling waves are not relevant here and 
have been neglected (i.e. k = 0) in writing the above solution that is still composed of two asymmetric nonlinear 
traveling waves. The nonlinear symmetry breaking depends on the set of parameters in which the DOPO operates 
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FIG. 1: Dependence of R 2 on the pump parameter Eq above threshold (Eq = 1). Numerical data obtained from integration of 
Model Q (continuous line) compared with predictions of Eq. JSJ (dashed-dotted line) and Eq. @ (dashed line). The insert 
shows the ratio between the continuous and the dashed-dotted lines. The parameters are: 7 P — 7 S = 1, A p = 0, A a — — 1, 
a p = 0.25, a s = 0.5, and a s = 0.25. 
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FIG. 2: Dependence of R 2 on convection (a s ). Numerical data (continuous line) are compared with Eq. © (dashed-dotted 
line) and Eq. @ (dashed line). E = 1.05-Eg (dark lines), E = l.LEg (light lines), J P = -fs = 1, A p = -0.2, A s = -0.5, 
a p = 0.5, a s = 1. 



via the ratio between the intensities (i.e. R 2 = \A 2 (k c )\/\A 2 (—k c )\ = I(k c ) / I(—k c )) of the two transverse modes of 
the signal [Eq. . This ratio has the explicit form 

with ui = 1 - [Re(F 3 )/2Re(/3)](/i - 1) and u 2 = [Im(F 3 )/2 Re(/3)]{fx - 1). 

This is the main analytical result. It allows a quantitative characterization of the nonlinear convection effects. 
Equation J5j emphasizes the coupling between convection and the distance from threshold. In the absence of convection 
R = 1, the two transverse modes have the same intensity and the amplitude equation exhibits stationary rolls. The 
presence of convection greatly complicates the expression of R 2 . However, near threshold (fj, > 1), the ratio of 
intensities R 2 , up to the leading order in [i — 1, is 

C(l + A 2 ) 

fl2 - 1 + ^en + (l + V + C2+ ^) ^- 1 ) ^ 

Note that R 2 — 1 is an odd function of a s , reflecting the importance of the sign of the velocity convection. Therefore, 
the choice of the convection direction (±a s ) can be useful to select one of the two modes by enhancing its parametric 



gain. Figure shows a typical variation of R 2 upon the physical pump amplitude Eq = My 1 + A 2 . We find a very 

good agreement between the analytical ratio R 2 and the numerical simulations of the Eq. In order to set the 
validity range of our predictions we have plotted the results obtained by increasing the pump till twice the threshold. 
Even for pump values 20% above threshold the agreement is within l%o(see insert in Fig.QJ. 

The ratio between the intensities of the two critical modes also provides the quantitative characterization of convec- 
tion in the nonlinear symmetry breaking. The numerical and analytical estimation of R 2 versus convection parameter 
a s , displayed in Fig. [21 & re again in very good agreement. Finally, note on Fig. |2|the existence of extrema leading to 
the most asymmetric configuration. 
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FIG. 3: Analytical velocity given by Eq. Q (continuous line) compared with numerical results (symbols) from the original 
model fl) . In the inserted plot we compare the velocities far from threshold. Same parameters of Fig. Q 



Let us now comment on the physics underlying the nonlinear symmetry breaking induced by convection. The most 
relevant physical parameter in the nonlinear interaction, above threshold, stems from the difference in frequencies 
of oscillations of each mode (±fc c ). This difference results from the presence of convection and disappears with it. 
Although, both traveling waves are propagating in the direction of the convection, their phase rotation arc no more 
opposite. Hence, the two traveling modes interact with a time delay with the pump. This gives rise to different gain 
from the pump leading to the nonlinear symmetry breaking observed in the signal. The energy transfer depends on two 
time scales and thus involves the pump decay rate ("f p ) as can be seen from the expression of D. We emphasize that in 
contrast with all previous studies dealing with the weakly nonlinear dynamics of OPO near threshold, 7 P appears, for 
the first time, in the cubic Ginzburg-Landau model because of the induced pump excitation phenomenon. This fixes 
the parameter range of the pump decay rates leading to a nonlinear symmetry breaking in the generated traveling 
waves. The stronger the pump decay rate, the weaker the asymmetry is. In the limit of adiabatic elimination of 
the pump, no asymmetry exists in the signal, consistently with the possibility to remove the walk-off by a change of 
reference frame. We have performed numerical simulations (not shown) with the same parameters as in Fig. ^ except 
that 7 P is decreased ten times. In this case we have observed a vanishing asymmetry (R 2 ~ 1) with respect to the 
result of Fig. ^ 

The second important feature that results from the convection induced nonlinear symmetry breaking concerns 
the propagation velocity of the generated traveling waves. Indeed, the convection effect on the signal is not only 
a translation of its transverse profile at the convection velocity. An increase in the pump enhances the action of 
convection in the nonlinear coupling between the fields leading to the velocity variation with the pump intensity. So 
that, if we set j s n = w cor , the corrected frequency of the traveling waves is ojr — uj c + u cor . Therefore their actual 
velocity is given by 

Uc + ls^ Im(/3) . 

- -a s j s -js -, „ rk O"- 1 ) ( 7 ) 



k c ,a fc c Re(/3) 

The above velocity expression shows that, in addition to the usual convection velocity (the first term of the right 
hand side) there is an excess velocity depending on the convection but, and most interestingly, it depends linearly 
on the incident pump above threshold (fi — 1). Figure [3] shows the predicted deviation of the actual velocity from 
the velocity convection by integrating the full nonlinear equations governing the DOPO dynamics [Eqs. ifTpt] when 
increasing the pump Eq till twice the threshold. As can be seen from this figure, there is a very good quantitative 
agreement for a pump till 10% above threshold. Only at threshold the nonlinear waves velocities coincide with the 
convection velocity. 

To summarize we have shown, in case of a degenerate optical parametric oscillator, that convection (walk-off) 
induced a nonlinear symmetry breaking in the traveling waves. We have also demonstrated that near threshold this 
mechanism is still described by a Ginzburg-Landau model with an original dependence of the nonlinear self-coupling 
term upon convection. As a result, nonlinear traveling waves are no more symmetrical and the explicit analytical 
expressions of their intensities variations with both convection and the distance from threshold are derived. Moreover, 
convection leads to nonlinear phase modulations that give rise to an interesting variation of the traveling waves velocity 
with the distance from threshold. Besides the context of optics, our results are relevant to many spatially nonlinear 
extended systems with convection. For instance, in the context of hydrodynamics, the competition between right- 
and left- propagating nonlinear waves in the convective flow, generated by a horizontal thermal gradient, leads to 
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an experimental observation of the nonlinear symmetry breaking. The broken symmetry has been evidenced via the 
estimation of the variation of the amplitude ratio of the right and left waves with the distance from threshold . 
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